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1. Introduction 

Photonic crystal fibres (PCFs), whose cladding is composed of a two-dimensional (2D) 
photonic crystal [TJ [2] , may confine and guide light through either a photonic bandgap 
effect jSl E] or by an effective high refractive index guiding mechanism Ej- Both 
classes of fibres have been the subject of numerous research and for a review we 
refer to Refs. [3 [H] and references therein. The latter class of PCFs has attractive 
features such as, broad-band single-mode (SM) operation [9 , possibilities for dispersion 
engineering PUJ [JT], and tailorable mode area ^2]- The 2D photonic crystal of the 
cladding not only provides more design freedom on engineering basic properties of the 
fibre, but also broadens the potential application of PCFs by the freedom to remove 
more air holes from the cladding f3J IE] or by introducing additional materials to 
the air holes [TBI 113 CHI EE ED] , both of which facilitate novel device operations based on 
PCFs. In the present work, the "PCF device" term is restricted to the ones obtained by 
removing some holes from the cladding such as, PCF couplers fSl E] , PCF polarization 
beam splitters [To] , or PCF distributed Bragg reflectors [2*T] . 

In the present paper, we describe an approach, which relies on modal transmission- 
line theory (MTLT), for three-dimensional (3D) investigations of the propagation of an 
optical beam launched into a PCF device. According to our knowledge only the finite 
element beam propagation method (FE-BPM) [22] has been adopted and utilized for 
doing such simulations. The FE-BPM is numerically robust, versatile, and applicable 
to a wide variety of structures. Unfortunately, this is often achieved at the expense of 
long computational times and large memory requirements, both of which can become 
critical issues especially when structures with large dimensions are considered or when 
used within an iterative design environment. 

MTLT, which has been developed for modelling multi-layered periodic media [23H2UI23 
ESlEZj- has been used for analysis of distributed feedback (DFB) lasers [23 EE], quantum 
well infrared photodetectors (QWIP) [2H], holographic power splitter /combiners [2H] 
and grating assisted integrated optics devices [HO]- Recently, it has also successfully 
been applied in a study of radiation fields from end-facet PCFs [HT]. MTLT relies on 
a plane-wave expansion of electromagnetic fields in the periodic media. Interpreting 
the plane waves as transmission-lines provides a systematic framework for study of 
wave propagation in multi-layered periodic media. Besides that, one can exploit all the 
concepts and methods of transmission-line theory [32] and electrical network theory [HS] 
for the study of wave propagation. MTLT has recently been developed for modal analysis 
of arbitrary shape optical waveguides |3lj . Here, we add a novel approach to this theory 
and utilize it for a three-dimensional study of propagation of light in photonic crystal 
fibre devices. 

The remaining part of the paper is organized as follows. In Section 2, we give a 
brief account of MTLT and describe the approach we use. In Section 3 we investigate 
examples that illustrate the utilization of this approach in the modelling of PCF devices. 
Finally, conclusions are given in Section 4. 
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2. Formalism 



The typical PCF device that we have in mind is composed of J layers with different 
relative permittivity functions e r j(x, y), j = 1, 2, ■ • • , J, which is illustrated in figure [T] 
We define Zj clS cl convenient local coordinate obeying < < tj and we 

consider wave propagation along the longitudinal direction of the structure, i.e. z-axis. 
Throughout this paper we consider non-magnetic materials with relative permeability 
fi r = 1 and all electromagnetic fields have a harmonic temporal dependence, exp(iujt). 
PCF devices with the typical shape shown in figure d usually confine light within the 
fibre core or cores. However, some application such as long-period fibre Bragg gratings 
do not share this feature so in the present study we exclude these applications. For the 
applications with spatially localized modes we use a super-cell approach and repeat the 
structure in the transverse xy-plane, along x and y directions with periodicities of T x 
and T y , respectively. It is assumed that the periodically repeated devices are separated 
by a sufficient amount of background region, here microstructured cladding, that their 
electromagnetic fields do not affect each other significantly. 

We want to study a PCF device when it is illuminated with an incident field 
Ei n c(r) propagating in layer 1 of the structure shown in figure Q This incident field 
is a solution to the source-free Maxwell equations in the PCF with the refractive index 
profile of layer 1. In this section we briefly address MTLT and modal analysis of optical 
waveguides using this theory. Subsequently we describe an approach, based on MTLT, 
for investigating the scattering and propagation of light in PCF devices. Throughout 
the paper, vectorial components are denoted by an arrow placed above them. The bold- 
style notation with uppercase and lowercase characters is used to designate matrices 
and vectors, respectively. 



2.1. Modal Transmission- Line Theory 

Embody a periodic medium with permittivity variation eoe r (x, y) = eoe r (x + T x ,y + T y ) 
and permeability The permittivity can then conveniently be expressed in the form 
of a two-dimensional Fourier series 

M n 2 ^ 2 ^ 

e r (x,y) = ^li^ Yl gm « ex P {- im Y x ) exp (~ in Y y ) 



where 

rT, 



— I [ e r {x,y)exp (im^-x) exp (in-^-y) dydx 

'^y Jo Jo ^ J-x ' ^ J-y ' 



T ^ 



(2) 



The electromagnetic fields must of course reflect the periodicity of e r (x, y) and according 
to the Floquet-Bloch theorem the fields in the doubly periodic medium are pseudo- 
periodic functions 

A(r) = lim V" V" a mn (2;)exp -if K x0 + m— )x exp -i( K y0 + n— )y (3) 

M,N — >oo '—^ *■ — ' L\ It / J I \ i,,/- 



m=-M n=-N 



Dabirian et al. : Propagation of Light in Photonic Crystal Fibre Devices 



4 



where, K = K x0 x + K y0 y is the Bloch wave-vector and A can be any of the 
electromagnetic fields E, H, or D. In order to facilitate calculations in matrix form, we 
introduce e, h, and d vectors whose elements are e mn , h mn , and d mn , respectively. The 
dimension of each vectorial component of the e, h, or d vectors in Cartesian coordinates 
(i.e. e x , e y , e z , h x , etc.) is 1 x (2iV+ 1)(2M + 1). Using these vectors, the constitutive 
relation D = eo€ r E converts into d = e Ne, where N is a square matrix whose elements 
are e mn and they are arranged in N in such a way that the equality d = eoNe holds. 
The temporal harmonic electromagnetic fields in a dielectric medium are solutions of 
the following source-free Maxwell equations 

V x E(r) = —itofXoH(r 

V x H(r) = iuD{f 

Using ((TJ), (jSJ), and vectors e, h and d in the source-free Maxwell equations these 
equations are transformed into the following system of differential equations: 
dv 



(4) 



dz 



di 

dz 



(5) 



#i = -iuCv 



or 



d\ 
~d? 

dH 



-cu 2 LCv 



(6) 



dz 2 



-u 2 CU 

where L and C are obtained in the calculations [2E! and 



(7) 



Equation (0) has the well-known form of telegraphist's equations for a multi-conductor 
transmission-line and we have emphasized the analogy by the choice of symbols 
so that e.g. i and v are interpreted as effective currents and voltages, respectively. 
Likewise, inductance and capacitance matrices of the multi-conductor transmission-line 
are denoted by L and C, respectively. In equation ©, u> 2 ~LC and u 2 CL are matrices 
with non-zero off-diagonal elements. We can formally diagonalize w 2 LC and uj 2 CL 
matrices using relations w 2 LC = PK 2 P 1 and w 2 CL = QK 2 Q _1 , where K 2 is a 
diagonal matrix. The diagonal elements of K 2 are eigenvalues of w 2 LC or w 2 CL. Here, 
P and Q are matrices whose columns are the eigenvectors of their relevant non-diagonal 
matrices. Once the K 2 and P have been determined, the matrix Q is also given by 
wCPK- 1 . 

From the above discussion it follows that (jOJ may be transformed into 
' d 2 w 



dz 1 



dH 



K 2 ^ 



(8) 
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where 

v = Pv , i = Qi. (9) 

In this new basis the transmission-lines are uncoupled and one may, in analogy with 
conductance eigen- channels in quantum transport [33], think of these new lines as the 
eigen-lines of the transmission-line system. Wave propagation in the periodic medium 
is described by K 2 , P, Q, see © and ©. 

Evidently from MTLT, the w 2 LC describes the propagation characteristics of 
longitudinal space harmonics. Eigenvalues of this matrix specify the square values of 
propagation constants of space harmonics. The propagation constants are obtained from 
the diagonal matrix K 2 considering the following condition |2lj : 

lm(K k ) + Re(K k ) < 0. (10) 

Electromagnetic fields of each space harmonic with a specified propagation constant are 
determined from its relevant eigenvector. 



2.2. Equivalent Network of Multi-Layered Media 

Consider the typical structure shown in figure The modelling task begins by 
periodically repeating the device in the transverse xy-plane with sufficiently large 
periodicities. As discussed above, wave propagation in each layer of this periodically 
repeated structure could be modelled by a transmission-line network whose behavior 
is described by (JHJ). Schematically, the equivalent transmission-line network of the j-th 
layer of this structure is depicted in figure El (a). In this figure the box containing Pj, Qj 
represents the consideration in (jOJ). 

A concise and effective formulation of voltages and currents of this transmission-line 
network can be described by: 

{vj = exp[-zKj(^ - ^)]vj,inc + exp[zKj(^ - tj)}^ 
ij = exp [-zKj (z, - tj )] ij, inc - exp [iKj {z s - ^)]ij,r 

where Vji nc , ij,inc, Vjr, and ij ir are vectors for the incident voltage, incident 
current, reflected voltage, and reflected current, respectively. The Kj = 
diag(ifji, Kj 2 , • • ■ , Kjk, ■ ■ •) is a diagonal matrix obtained by computing the square root 
of the K 2 matrix. The exp[— zKj(tj — Zj)} is also a diagonal matrix with diagonal 
elements exp[— Kj k (tj — Zj)]. 

Essential electromagnetic boundary conditions could be simply satisfied at the interface 
of two different layers by the continuity of voltages and currents in transmission-line 
theory. At the interface of typical different I and / + 1 layers, illustrated in figure El (b), 
the continuity rule is described by: 

p j^tj = Pj+iyoj+i ^ 
Qjitj = Qj+iioj+i 
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where v t j, v v t j, and v j +1 have been defined in figure|21 On the basis of MTLT the 
transmission-line network of the periodically repeated typical PCF device is illustrated 
in figure El In the equivalent network of figure El and also in numerical simulation a 
total height of hj is considered. At the beginning (z = 0) and the end (z = hj) of 
the structure, the well-known radiation condition of electromagnetic theory is applied, 
which is depicted in the equivalent network by match load. Here we exploit a primary 
feature of radiation condition; i.e. the zero reflection at these points. 

2.3. The Approach 

Consider the structure shown in figure ^ When this structure is illuminated with an 
incident electromagnetic field Ei nc (r), propagating in layer 1 along the positive direction 
of z-axis, the total field in layer 1, Ei(r), is given by 

£x(f) = E inc {r) + E^ r {r) (13) 

where Ei jr (r) is the reflected field inside layer 1. The incident field is usually a fibre 
mode so for investigating its interaction with other layers we must calculate it and then 
calculate Ei >r (f) and finally E\{f). From the known field at layer 1, we calculate the 
fields in other layers utilizing equations (fTTj) and (pH2|) . 

Since the incident field E inc (f) is a guided mode of the waveguide with refractive 
index profile of layer 1, we could determine it utilizing MTLT, exploiting the features 
of transmission-lines This calculation is achieved by examining the out of 

plane propagation of a periodic medium whose refractive index variation is obtained 
by periodically repeating the waveguide in transverse plane with sufficiently large 
periodicities. Evidently from transmission-line theory, the matrix a; 2 LC contains the 
information of the out of plane propagating waves, called space harmonics in the field 
of diffraction grating [23J. Eigenvalues of the matrix w 2 LC, diagonal elements of K 2 , 
specify squared-values of propagation constants of these space harmonics. Each column 
of the matrix P, an eigenvector of the matrix w 2 LC, describe the electric field profile 
of its relevant eigenvalue. Among the space harmonics the ones whose field profiles are 
localized within the waveguide specify guided modes of the waveguide. In index guiding 
waveguides this condition is simplified to the refractive index guiding condition. 
Through the modal analysis of the fibre with layer 1 refractive index profile we determine 
Vj,inc(-2 = 0). Afterwards, for complete determination of field in the first layer, 
calculation of Vj ir (z = 0) and ij, r (z = 0) is also required. These values are obtained 
by the following relations: 

v j>r (z = 0) = R&vj.U* = 0) 

ij, r (*=0) = RJ^incO = 0) 

where R^ is the upward reflectance matrix at z = 0. Generally we define R"j as the 
reflectance matrix of a propagating wave along the positive direction of z-axis at the 
local geometry Zj] for instance Rq X is the upward reflectance matrix at Zj = for j = 1. 
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The variation of R"j along z is treated by the following relation [H7j 

{Rtj = F r j(R j+i) 
R j = exp(Kjx)R t j exp(Kjx) 

where F rj (R j +1 ) is composed of the following set of equations: 



(15) 




'tj - 



tj - 




(Z tj -I)(Z tj + I)- 



(16) 



Computation of Rq X is started from the topmost layer, where the reflectance is zero. 
Considering (JTHJ) and fT6j) at each layer and interface of layers, the Rqx would be 
calculated. From the known Rq 15 the vi r and ii ir are determined by Equation (JHJ). 
Electromagnetic fields at other points of the first layer would be computed using (JT3J). 
Inside other layers, electromagnetic fields will be calculated using (fTT[) and (fT2"|) . 

3. Validation and Numerical Implementations 

In this section, several examples will be considered to illustrate and also validate the 
utilization of the proposed approach. 

3.1. PCF coupler 

The cross-section of the PCF-coupler we want to study is depicted in the inset of figure^] 
It is composed of a triangular lattice of air-holes in silica with two missing air holes. We 
validate the described approach in the present paper by verifying the coupling action 
of the coupler and comparing the obtained coupling length through this approach with 
the one obtained by considering even and odd modes. In the simulation the pitch, 
A = 7.2 /im, and normalized hole-diameter to the pitch, d/A = 0.45 have been set. We 
perform the simulation at the normalized wavelength A/A = 2txc/(ujA) = 0.1. 
In simulation it is assumed that the light is launched into one core of the coupler, for 
instance core A, by butt-coupling of a similar single-core fibre whose core is aligned to 
the core A. The coupler and the fibre coupled to it constitute a two layer medium, which 
could be considered an example of the general case of figure [TJ 

As it is described in Section 2.2, at first we repeat the structure periodically in 
the transverse xy-pl&ne with 10A x 10A periodicity. Fiber cores of both the single-core, 
first layer, and the double-core, second layer, are considered as defects so treated by 
the supercell approach [36J. We calculate the fundamental mode of the single-core fiber 
using a MTLT-based approach of [31] which has been briefly described in Section 2.3. 
Through the simulation we obtain the fiber mode as the voltage and current vectors 
Vj,inc(<2 = 0) and ij ( in C (^ = 0). These vectors describe the electromagnetic fields of the 
fiber and are related to the fields through the Eqs. (7) and 0. Evidently from (J3J 
the fiber mode is the weighted summation of individual plane waves with different wave 
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vectors. From the known vj ^ nc (z = 0) and ij t i nc { z — 0) electromagnetic fields inside all 
the structure will be computed by tracking the approach described in Section 2.3. 

We illustrate in figure El the normalized electric field intensity when the HEu mode 
of the single-core fibre, travelling across the z-axis, is launched into the core A of the 
dual-core fibre. Inside the dual-core fibre, the light starts coupling from the core A to 
core B. Up to the distance of 1440 /zm from the interface of the coupler and single-core 
fibre (z = 2.44 cm) all the confined light in the core A will be coupled to the core B. 
This distance is called coupling length, L c , and alternatively may be computed from 
the difference of the propagation constants of even, /3 e , and odd, (3 a , modes of the dual- 
core fibre through the relation of L c — ir/ | /3 e — (5 Q |. The computed coupling length 
between the even and odd modes in translational invariant system is 1410 /im, which 
is in a good agreement with that obtained through the approach of this paper. The 
normalized intensity of electric field in the center of core A is depicted in figure EJ The 
coupling length has been indicated on the figure. 

3.2. PCF Bragg grating 

The case of a PCF Bragg grating arises in various advanced applications of photonic 
crystals. In PCF lasers an optical cavity may be formed through two PCF Bragg gratings 
created by introducing a spatial periodic modulation of the refractive index to the fibre 
core along the fibre axis [21]. Photonic crystal vertical cavity surface emitting laser (PC- 
VCSEL) |3H] is a novel application of the photonic crystal in the laser application, which 
is similar to standard VCSELs except that a photonic crystal structure is defined by 
introducing regular lattice of air-holes with one missing air-hole to the top mirror. These 
lasers, as well as the single mode operation, have side-mode suppression ratios about 
35-40dB jsHJ. These attractive features are facilitated by the presence of the regular 
lattice of air holes as has been studied qualitatively utilizing concepts of PCFs [39J. 
Using the approach described in this paper, the reflectivity from the top mirror could 
be investigated three-dimensionally. The modelling of the laser mirrors is generally a 
crucial issue in the design and analysis of lasers jlU]. Here we analyze an example of 
an in-plane grating in a PCF to illustrate the proposed approach. The structure under 
consideration is depicted in figure EH The cross section of each layer is a square-lattice 
photonic crystal composed of air holes in the background material with one missing 
air hole. The layers specified by white color have refractive index 1.45 and the colored 
ones have refractive index 1.6. The air-holes of normalized diameter d/A = 0.53 are 
arranged on a square lattice with pitch A = 7.2 /mi. Such mirrors have recently been 
utilized as the top distributed Bragg reflector of PC- VCSELs |JT]. In the structure 
the thickness of colored layers, t, is 0.12/zm and the periodicity of the Bragg mirror 
is a = 0.245 /im. Utilizing the approach of this paper, we examine the interaction 
of the travelling fundamental mode of the first layer with the grating at A/A = 0.1. 
Figure shows the two-dimensional intensity plot of the electric field in the case where 
the fundamental mode of the squared- lattice PCF (with lattice index of 1.45) is incident 
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on the mirror. The incident field is partially reflected at interfaces of different layers, 
leading to an interference pattern caused by interference of the incident field and the 
reflected ones. Figure [3 also illustrates how perfectly boundary conditions at different 
material interfaces of the distributed Bragg mirror are fulfilled. 

4. Conclusions 

Optical properties of PCFs may typically be successfully analyzed within the assumption 
of translational invariance along the fibre axis. However, in real life the important device 
applications employ PCFs of finite length and the hypothesis of translational invariance 
is not applicable. In this work we have described a semi-analytical approach for three- 
dimensional fully vectorial analysis of photonic crystal fibre devices. Our approach 
rest on the foundation of modal transmission-line theory and offers a computationally 
competitive alternative to beam propagation methods. The approach is illustrated by 
simulations of the coupling action in a photonic crystal fibre coupler and the modal 
reflectivity in a photonic crystal fibre distributed Bragg reflector. 
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Figure 2. Equivalent electrical networks elements (a) transmission-line unit presenting 
a single layer (b) General junction of two transmission-line units at different layer 
interface. 
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Figure 3. Equivalent transmission- line network of the multi-layered structure shown 
in Fig. HJ 




Figure 4. Intensity of electric field in the center of the core A. The inset shows the 
cross section of the PCF coupler with the two cores A and B. 




Figure 5. Distribution of electric field intensity (a) in the xz plane at the centers of 
the fibre cores (b) in the xy plane. The electric field in the xy plane is computed at 
the center of the core A. 




Figure 6. Geometry of the photonic crystal fibre distributed Bragg reflector. 




Figure 7. Distribution of electric field intensity (a) on the xz plane at the center of 
the fibre core and (b) at the center of the fibre core. 



